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Abstract. A new type of quantum transfer matrix, arising as a Cholesky factor for 
the steady state density matrix of a dissipative Markovian process associated with 
the boundary-driven Lindblad equation for the isotropic spin- 1/2 Heisenberg (XXX) 
chain, is presented. The transfer matrix forms a commuting family of non-Hermitian 
operators depending on the spectral parameter which is essentially the strength of 
■^j- ■ dissipative coupling at the boundaries. The intertwining of the corresponding Lax and 

ON . monodromy matrices is performed by an infinitely dimensional Yang-Baxter i?-matrix 

which we construct explicitly and which is essentially different from the standard 



4x4 XXX i?-matrix. We also discuss a possibility to construct Bethe Ansatz for 
the spectrum and eigenstates of the non-equilibrium steady state density operator. 



Furthermore, we indicate the existence of a deformed i?-matrix in the infinitely- 
dimensional auxiliary space for the anisotropic XXZ spin-1/2 chain which in general 
provides a sequence of new, possibly quasi-local, conserved quantities of the bulk XXZ 
dynamics. 
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1. Introduction 

The theory of integrable quantum systems in 1 + 1 dimensions, the so-called quantum 
inverse scattering, is a well developed field of mathematical physics [30J [32j |9] [16] which 
pioneered important new algebraic structures in pure mathematics, such as quantum 
groups and their representations. The fundamental object in this theory is the i?-matrix, 
a solution of the celebrated Yang-Baxter equation [2], which gives rise to integrable 
Hamiltonians possessing infinite families of conserved quantities. Furthermore, these 
techniques often lead to explicit methods for diagonalizing the Hamiltonian, such as 
Algebraic Bethe ansatz (ABA) [30| |9j [16] or Baxter Q-operator [3]. More recently, 
the theory of integrable quantum systems also found applications in classical non- 
equilibrium physics, namely in solving markovian stochastic many-body interacting 
systems such as asymmetric simple exclusion process [271 IS]- There has been even 
an attempt to develop a non-equilibrium Bethe ansatz approach to quantum impurity 
problems [19] , nevertheless the practical feasibility of these technique and its relation to 
general integrability structures such as Yang-Baxter equations remains unclear. 

However, very recently explicit results appeared for driven quantum many- 
body systems with a strong interaction, namely a closed matrix product ansatz 
(MPA) for non-equilibrium steady state (NESS) density-operator of the boundary- 
driven Lindblad equation [221 E31 El] of an anisotropic Heisenberg (XXZ) spin 
1/2 chain. Lindblad equation [HI [11] is the canonical model of continuous-time 
Markovian quantum dynamics. This solution has been later interpreted in terms 
of infinite-dimensional representations of Lie algebra sl(2), and its quantum-group 
deformation for the anisotropic spin interaction, and generalized to more general 
boundary dissipators/drivings [15]. Remarkably, perturbative expansion of NESS in the 
dissipation strength gave rise to a novel XXZ quasi-local conservation law [22] which is 
unrelated to previously known local conserved quantities of the XXZ chain [TO] derived 
from the 'standard' XXZ transfer- matrix, and which has important consequences for 
understanding ballistic transport at high-temperatures [13] . 

In this paper we put these results [221 [231 [2U [TO] into the framework of the theory of 
integrable systems. Focusing mainly on the isotropic case (XXX model), we rigorously 
construct an R— matrix satisfying Yang-Baxter in an infinitely dimensional auxiliary 
space which carries irreducible infinitely-dimensional representation of sl(2), so that the 
corresponding family of commuting transfer-matrices is given by the Cholesky factor 
of the unnormalized NESS density operator [23]. However, the commuting transfer 
matrix is given as the ground-state matrix element of the monodromy matrix, and 
not as its trace as in standard ABA, and is neither a Hermitian nor a diagonalizable 
operator, which is a manifestation of far-from-equilibrium character of the problem. As 
the spectral parameter in our R— matrix comes from the boundary dissipative coupling 
we chose to call our formalism the exterior integrability. One may also provide arguments 
for the existence of a deformed version of the infinite dimensional exterior R— matrix 
in the anisotropic (XXZ) case. Two important immediate applications of the new 
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R— matrix are proposed: (i) Construction of an infinite family of new almost-conserved 
[13] quantities mutually in involution which include the one discussed in [22] and which 
should shed further light on the understanding of finite-temperature quantum transport 
problem [291 [23 H2], and (ii) Construction of ABA for diagonalization of NESS density 
operator. We stress that even if the exterior integrability is defined with respect to 
particular integrable dissipative boundaries, it may produce interesting new results 
for bulk properties of the system in the thermodynamic limit, such as the quasi- local 
conserved quantities. 

After defining the main concepts of non-equilibrium quantum integrability of 
the XXX model in Section 2, we write explicit expression for the corresponding 
infinite R— matrix in Section 3 and rigorously prove that it satisfies the Yang-Baxter 
equation. In Section 4 we describe some interesting properties of the i?-matrix and the 
corresponding non-equilibrium monodromy matrix. In Section 5 we discuss potential 
applications and extension to an anisotropic case, and conclude. Some technical aspects 
of our proofs are put into appendices. While material presented in Sections 2, 3 and 
the appendices should be mathematically rigorous, further results discussed in Sections 
4 and 5 are partly based on heuristic and empirical arguments. 

2. Exterior integrability of the nonequilibrium steady state 

We focus on the stationary Lindblad equation for the NESS density operator p^ 

i[Hxxx,Poo} = eV( Poo ) (2.1) 

for the Heisenberg XXX Hamiltonain of a chain of n spins 1 /2 

n-X 

H X xx = y^ l2J-i®fa8l2»-:*-i, h = 2o- + <g>cr~+2o-~<g>o- + +o- z <g>cr z , (2.2) 
i=i 

where o~ ± = |(o" x ± ic" y ), o" z , c = I2 are standard Pauli matrices acting over a 

2— dimensional quantum spin space R s ~ C 2 and 1^ is a <i-dimensional unit matrix. 

We chose the simplest solvable far-from-equilibrium dissipative driving [23J with a pair 

of Lindblad jump operators with dissipation-driving strength e 

2 

V(p) = J2( 2L kP L l-{ L l L k,p}), L 1 = a + ^l 2 n-i, L 2 = l 2 n-i®(j-.(2.3) 
k=i 

As it has been shown in [23], the unique NESS density operator can be written explicitly 
in the Cholesky factorized form 

Poo = S(X)S\X), Poo = -^L. (2.4) 

trpoo 

where the operator S(X) admits an elegant representation in terms of MPA: 

S(A) = J2 (0|A Sl (A)---A Sn (A)|0)(7 Sl ®---®a s " (2.5) 

si,...,s n e{+,o,— } 
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where 

oo 

A (A) = j>2(A)|*)<*|, 

fc=0 
oo 

A + (A)=5>+(A)|fc><fc + l|, (2.6) 

fc=0 
oo 

A-(A)=^%(A)|A; + 1)(A;|, 

fc=0 

is family of tridiagonal matrix operators acting on an infinitely- dimensional auxiliary 
Hilbert space H a with a canonical basis {|0) , |1) , |2) , . . .}. In fact, the consistency of 
solution (I2.5P with the defining equation (j2.ip requires that the matrix operators span 
an infinitely dimensional irreducible representation of sl(2) algebra 

[A + (A), A_(A)] = -2A (A), [A (A), A ± (A)] = ±A±(A) (2.7) 

which may be - up to unitary transformations - uniquely chosen asjl 

a° k (\) = \-k, a+(A) = £;-2A, a -(\) = k + l, A G C, (2.8) 

with the complex representation parameter A being fixed by the boundary dissipation 
strength 

A = -. (2.9) 

Defining a A-dependent linear operator from End("H s <g> "H a ) 

L(A) = <x°®A (A)+a + ®A + (A)+a~®A_(A) = ( 

the Cholesky factor can be expressed even more elegantly [15] 

S(X) = (0| L(A) ® s L(A) ® B . . . ® s L(A) |0) = (0| L(Af sn |0) . (2.11) 



A (A) A + (A) 
A_(A) A (A) 



(2.10) 



Here and below we use the following compact and unambiguous notational convention. 
For operator-valued matrices, we use a symbol ® s as a partial tensor product, namely 
it implies a tensor product with respect to the quantum spin space % s and an ordinary 
operator/matrix product with respect to the auxiliary space % & . Analogously, £g> a will 
denote a tensor product with respect to "H a , and a matrix product in 7-L s . For example, 
for X G End(?C <8> U® k ), Y G End(Hf ® U® m ), X ® s Y G End(^ s ®°' +0 <g> -H® fc ) 
making sense if k = m, and X ® a Y G End("H® : ' ® ?^ a ) making sense if j = I. To 

emphasize the exterior integrability concepts we shall write in bold all symbols which 
are not scalars with respect to auxiliary space % a . 

The key step of this work is to recognize that L(A) can be interpreted as the Lax 
matrix (the so-called L— matrix) matrix of a novel integrable system. This is founded on 

| Note that the representation of Ref.[53], specialized to the isotropic case, is, up to a constant, unitarily 
equivalent to ([2] 
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a simple empirical observation, namely that the Cholesky factors commute for arbitrary 
complex values of the representation/dissipation parameters 

[S(X),S(fi)] = 0, VA,/iGC. (2.12) 

This observation can be understood as a consequence of existence of an R— matrbq§| 
R(A, /x) G End("H a <8> "H a ) for almost any A,/i G C, to be shown in Section 3, which 
satisfies the so-called RLL (or local intertwining) relation: 

R(A, fi) (L(A) ® a L(//)) = (LGu) ® a L(A)) R(A, /i), (2.13) 

Following the procedure of ABA [16] the local intertwining relation immediately implies 
intertwining for a product of the so-called monodromy matrices T(A) G End( - Hf n ® - H a ): 

T(A) = L(A) ® s L(A) ® B . . . ® B L(A) = L(Af s ™, (2.14) 

namely 

R(A, fi) (T(A) ® a T(/i)) = (T(jjl) ® a T(A)) R(A, //). (2.15) 

Indeed, Eq. fl2. 13|) implies Eq. (" 12. 15[) after noticing that, due to associativity of matrix 
multiplication: 

T(A) ® a T(/i) = (L(A)®- n ) ® a (L(//)®" n ) = (L(A) ® a L^))®-" (2.16) 

Unlike in the standard formalism of ABA where the auxiliary space is finite 
dimensional and the concept of a transfer matrix is usually associated to the partial 
trace of monodromy matrix with respect to the auxiliary space, we define here the 
auxiliary ground state expectation (0| T(A) |0) = 5(A) as the transfer matrix. In order 
to establish the commutativity of the transfer matrix we also require, besides the RTT 
relations f !2.15j) . the R- matrix to satisfy additional boundary conditions 

(0,0|R(A,/i) = (0,0|, R(A,/i)|0,0) = |0,0), (2.17) 

where \k,l) := \k)®\l) , (k,l\ := (fc|(g)(Z|. Eq. (j2.12p then follows straightforwardly, after 
writing the transfer-matrix product in 7/ a <%> 7i a , 5(A)5(/x) = (0, 0| T(A) ® a T(//) |0, 0): 

S(\)S(fi) = (0, 0| R(A, /i) (T(A) ® a T(/i)) |0, 0) 

= (0,0| (T(/x) ® a T(A)R(A,/i)) |0,0) = S(ji)S(X). (2.18) 

Despite the boundary condition (12.171) may seem a-priori unjustified at the moment, we 
shall show further on, that such a property naturally follows from the so-called ice-rule 
property of the i?-matrix. 

It is perhaps remarkable that the transfer matrix of our problem 5(A) is non- 
Hermitian, non-normal, and even non-diagonalizable operator. Using the MPA form 
(12. lip we can write its matrix elements in the spin basis {\u) = \vi, . . . , u n ) ; Vj G {0, 1}} 
of nf n , a z \u) = (-l) u \u), as 

(E'\ S(X) \v) = (0| A„ 1 _ K (A)A, 2 _^(A) • • • A^_„ n (A) |0> (2.19) 

§ We follow nomenclature of Ref. |16) here. Often in literature the term i?-matrix is reserved for an 
operator PR where P is permutation operator which swaps the auxiliary spaces. 
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so that the rules (0| A = A (0| and (0| A_ = 0, following from representation (12. 6p . 
imply the matrix of 5(A) to be upper triangular, 

n n 

(e'\ S(X) \k) = 0, if J2 v 'j T ' J > J2 V 3 T ~^ ( 2 - 2 °) 

3=1 3=1 

and having a constant diagonal 

(k\ 5(A) \u) = A n . (2.21) 

Consequently, all eigenvalues of 5(A) are equal to A™, but since 5(A) is not a multiple 
of the identity operator it must have a non-trivial Jordan decomposition, i.e. it must 
be non-diagonalizable. 

Similarly, we can write the quantum space matrix elements of the general 
monodromy matrix elements 

If (A) := (k>\ T(A) |*> (2.22) 

following the expression f |2.14j) in terms of MPA 

(A T k k '(\) \k) = (k'\ A U1 _ K (\)A U2 _ U ,(\) ■ ■ ■ A Vn _ K (\) \k) . (2.23) 

Tridiagonality of operators ( 12. 6p immediately implies a magnetization selection rule, 
namely (I2.23J) vanishes unless 

n 

k'-k + J2 u 3~ v 'j = °- ( 2 - 24 ) 

3=1 

This in turn implies that T^ (A) changes the z— component of magnetization by 2(k' — k), 

[M,Tl:\\)} = 2(k'-k)T k k '(\), (2.25) 

writing magnetization operator as M := Y^=i I2J- 1 ® °" z ® 12»-j- 

S.i. Ice-rule - i/ie particle conservation law 
Let us write out the i?-matrix in components 

CXD OO 

R(A,/i)= J2 Y,B^{X,n)\k } kri(lJ\. (2.26) 

k,k'=0l,l'=0 

We will show in the following section that the exterior R— matrix of the XXX model 
(and also for a more general XXZ model, see subsection 15. 2p obeys a selection rule, 
namely Rffi (A,/x) 7^ only if k + k' — I + 1'. This can be interpreted as a particular 
particle conservation (global U(l)) symmetry of the R— matrix, meaning that it should 
commute with the particle number operator 

[R(A,/i),N]=0, (2.27) 

00 

N = -(A (0) <g> 1 + 1 <g> A (0)) = alo+i. (2.28) 

Consequently, one can interpret the i?-matrix as a particle-number conserving scattering 
matrix of a system of auxiliary quasi-particles. Decomposition (12. 28ft suggests a natural 
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splitting of a tensor product of two copies of auxiliary space into a direct sum of 
eigenspaces of N 

oo 

•H a ®-Ha = 0^i Q) . (2.29) 

a=0 

As we see, there are a + 1 states \k,a — k) within each sector 7ia ■ Therefore, for 
any X G End(% a ® "H a ) which commutes with N, [X, N] = 0, we shall denote with 
upper-bracketed index a an (a + 1) x (a + l)-matrix component of decomposition 
X = ©^Lq-^ ■ F° r exam pl e > we shall often write the R— matrix in the so-called 
ice-rule form 

oo a oo 

R(A,//) = Y, E tfi&ti M" k ) &«" *l = 0R (a) (A,/i). (2.30) 

a=0 k,l=0 a=0 

As elements of Hi are scalars, any i?— matrix satisfying the ice-rule (I2.30p should 
trivially obey the boundary condition (j2.17p . 

3. Exterior R- matrix 

Here we shall write out and prove our main result, an explicit form of the infinitely 
dimensional exterior R— matrix which satisfies the defining RLL relations ( 12.131) . 

Theorem. A solution of the RLL (12.131) relation for Lax operator (j2.10p reads 

R(x+ \y,x- \y) = exp(yH(a;)), (3.1) 

for any x G C\|Z + , y G C. The generator H(x) admits a block decomposition according 
to the ice-rule, 

H(*) = 0H(*)(x), (3.2) 

O: 

with explicit form of the matrix elements 

<?m = tf^Q B-U- (*;'_" 1 )m*), * > « + i, 

1 a— fe— 1 

#$(*) = -2 E w*)> 2A; ^ ( 3 - 3 ) 

m=fc 

where we introduced simple-pole functions f m (x) := (x — m/2)~ 1 . 

Proof. We start by using ( 13. ip as an ansatz and reparametrize the RLL relation (12. 13ft 
in the center- of -mass and displacement spectral parameters, 

x=(\ + fi)/2, y = \-n, (3.4) 

namely 

exp(j/H(x))(L(x+ij/)® a L(x-|j/)) = (L (x-|y) ® a L (x+|j/)) exp(yH(x)), (3.5) 
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yielding the form where non-trivial dependence enters through the generator H(x), in a 
way which resembles a Lie group structure. Furthermore, we employ the fact that the 
Lax matrix L(x) has a simple linear dependence on the spectral parameter 

L(A) = L + AL', L := L(0), V : = -j-L(ar), (3.6) 

whence 

A(x, y) := L (x + \y) <g) a L (x - \y) 

2 

= L{x) ® a L{x) - | (L(x) ® a 1/ - L' ® a L(x)) - V - (1/ ® a L') 

=: A (x) - |A X - ^A 2 . (3.7) 

At this point we emphasize that the whole x-dependence is absorbed into zero-th degree 
component Ao(x), whereas Ai^ are matrices with constant (x independent) elements. 
In particular, Ai = L(x) ® L' — L' <g> L(x) = L <8> L' — L' <g> L . Writing the Weyl basis 
of End("H s ) as E u,v = \u) (u'\ and expressing 

V = (E 00 + E u ) ® 1 - E 01 ® B, 

B:=-^A + (x) = 2j>>(fc + l|, (3.8) 

we can write the three orders A 0) i,2 (\S.7h as operators over % s ® %& ® % a factoring out 
the components in the physical space 

A (x) = E 00 g> (Ao(x) ® Ao(x) + A+(z) g> A_) 
+ £ n ® (A (x) <8> Ao(x) + A ® A + (x)) 
+ £ 01 g> (A (z) g> A+(x) + A+(x) g> Ao(x)) 
+ £ 10 <g> (A (x) ® A + A <g> Ao(x)) , (3.9) 

Ax = E 00 ® (A (0) ® 1 - 1 <g> A (0) + B <g> A_) 
+ £ n (8) (A (0) <g> 1 - 1 ® A (0) - A ® B) 
+ £ 01 ® (A+(0) g> 1 - 1 ® A+(0) + B ® A (0) - A (0) <g> B) 
+ £ 10 ® (A_ (8)1 -1® A_), (3.10) 

A 2 = (£ 00 + £ n ) <8> (1 ® 1) -E 01 (8)(B(8)1 + 1®B). (3.11) 

After inserting proposed ansatz for the solution (I3.ip . we shall expand (13. 5p in 
terms of nested commutators - (i) we multiply (I3.5J) by the operator exp (— |H(a?)) 
from the left and from the right, and (ii) we utilize the defining Lie-group identity 
exp( &dx)Y = e x Ye~ x , where adx := [X, •], which brings ( 13. 5ft to an equivalent form 

exp (ij/ad H (*))A(x,j/) - exp (-|yad H (x))A(x, -y) = 0, (3.12) 

or 

sinh (lyadn^)) (Aq(x) - jy 2 A 2 ) - |ycosh (|y adn^)) A x = 0. (3.13) 

Expanding the hyperbolic functions we obtain a power series in y, which always exists 
in terms of finite matrix exponentials due to decomposition ( 13. 2p . Clearly, since the 
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expression above is an odd function in y, we find only odd orders nonvanishing. In the 
first order in y we have 

ad H (x)A (:r) = A 1; (3.14) 

while for an arbitrary odd order y 2l+1 with I > 1: 

ad HS Ao ( x ) " ( 2/ + X ) ad H(x) A i - 2/(2/ + 1) ad^ 1 ) A 2 = 0. (3.15) 

The relation (13. 14ft is in fact an infinitesimal RLL relation for y — > dy and will be in 
the following referred to as HLL relation. 

Next we show that an infinite sequence of operator equations (13.151) can be in fact 
reduced to only two equations. The first one is just the third order [ (13. 15ft for / = 1] 
after substituting adjure) A (:r) from HLL relation f)3.14[) : 

adH (a; )Ai = -3ad H ( :c )A2, (3.16) 

Then we subsequently use fj3.14j) and A3. 16|) to eliminate A (x) and Ai from the sequence 
(I3.15P for any / > 1, arriving at ad^ I 7 1 NA 2 = 0, for which a sufficient condition is 



ad^ } A 2 = 0. (3.17) 

We have thus shown that three simple y— independent equations, namely (I3.14p . 
and a pair (I3.16f3.17p imply validity of Eq. (" 13. 15[) for any /, and consequently of the 
full RLL relation for any pair of spectral parameters x, y for which H(x) exists, i.e. 
x eC\|Z+, yeC. 

The remainder of the proof is thus to verify identities (I3.14|3.16|3.17p which we 
formulate in two lemmas below. 

3.1. The HLL relation 

Lemma 1. The generator of the i?-matrix H(x) (13 .3p solves the infinitesimal RLL 
relation f l3TT4"P : 

[H(a;),Ao(x)] = A 1 , (3.18) 

for any x G C \ |Z + for which it is defined. 
Proof. Let us define a permutation map - homomorphism - over End(H a (g> 7i a ) 

7r a (X) = PXP" 1 = PXP, (3.19) 

where P is a permutation operator over "H a <g> "H a , acting as 

P\k,l) = \l,k), k,leZ + . (3.20) 

Permutation operator conserves the number of auxiliary excitations, hence it satisfies 
the ice rule 

P = 0pM Pff = 5 k+l , a . (3.21) 
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We may write shortly 7r a (a®b) = b®a. Then we define another map tt s over operators 
in the quantum spin space End("H s ), by 

n s (E uu ') = E 1 -"'' 1 -", (3.22) 

or equivalently vr s (cr°) = <r°, ii^a^ = a ± , vr s (cr z ) = — cr z , so the full parity map over 
H s ® "H a ® "H a is defined as 

7T = 7T S (g) 7T a . (3.23) 

It is important to note that the operators Aq^^ and the generator H(ir) are 
eigenoperators of the parity map, i.e. they have well defined parities [see ( 13.9113. 10113. lip ]: 

7r(A fe ) = (-l) k \ k , k = 0,1,2, 7T a (H) = -H. (3.24) 

Notice that H(x) operates trivially (i.e. as a scalar) in the physical space H s . The whole 
expression ( 13.181) is then an eigenoperator of it with eigenvalue — 1, 

[H(z), Ao(a;)] - Ai + 7T QH(x), A (x)] - Ai) = 0. (3.25) 

Let us now introduce the components in the quantum spin space, either in Weyl or 
Pauli basis, A s k , A™ G End("H a <8> 'Ha), namely 

1 

v,i/'=o se{o,+,-,z} 

The identity (13.181) to be proven then writes 
1 
Y, E™' ® ([H(x), Af (x)} - Af) = 0. (3.27) 

u,u'=0 

whereas the symmetry relation (13.251) . noting (13.241) . can be rewritten as 

(E m -E n )®(([H(x),A°\x)]-A™) - ([H(x),Al\x)]-A{ 1 )) = 0. (3.28) 

This means that out of four equations in End("H a <8> "Ha), implied by ( 13. 2 7ft . only three 
are independent, say the components 00, 01 = + and 10 = — . 

Furthermore, we apply the a— decomposition of the A& operators 
00 
AZ = 0A2* (3.29) 

where A k a '°' z = A k ° ± A k a ' n G End (Ha ) are (a + 1) x (a + 1) matrices, while 
A[ a)+ G Lm(ni a) ,Hi a+1) ), and A { k a) ~ G Un(Hi a+1) ,H a a) ) are (a + 1) x (a + 2) and 
(a + 2) x (a + 1) matrices, respectively. With a bit of patience one can derive explicit 
expressions from Eqs. ( 13.9113. 10113. Ill) , for the constant operators Ai and A 2 . Using a 
compact notation for a canonical basis of H& , \k) = \k, a — k), the only non- vanishing 
blocks/components are 



A (a)o = ^ 2 (a-2ife)|A;)(A;| 



fc=0 
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a-l 

+ J2 { 2 ( a ~ k ) \ k ) ( k + M ~ 2 ( k + 1) |fc + 1) ( k \ } (3-30) 

fc=0 
a a— 1 

A (a) z = ^ 2 (a-fc)|fc)(fc + l| + ^2(fc + l)|fc + l)(fc|, (3.31) 

fc=0 fc=0 

A (a)+ = J2 |(3fc - a) |Jfe) (]fe| + (3ife - 2a) |]fe) (jfe + 1| }, (3.32) 

fc=0 
a 

A (a)- =J2{(k + l)\k + l)(k\ + (k-a-l) \k) (k\ }, (3.33) 

fc=0 
a 

^{-2|fc)(fc|-2|fc)(fc + l|}. (3.34) 



fc=0 
a 

A {a)+ 



fc=0 

The full set of finite matrix equations which remain to be verified then reads: 

[H^(x),A^ )00 (x)} = A ( f )00 , (3.35) 

H {Q \x)A^ + (x) - A^ )+ (x)H^ +1 \x) = A[ a)+ , (3.36) 

U^ +1 \x)A^ ] -(x) - A { a) ~(x)H {a) (x) = A[ a) ~. (3.37) 

For this one needs to show that for all equations residua at the possible poles, x = p/2, 
p = 0, ...,a+ 1, match as well as the remainders. This is done in full detail in 
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3.2. Master- symmetry of the H -matrix 

Lemma 2. For any x G C \ |Z + for which H(s) (I3.3P is defined, it satisfies identities 
(l3~TBT3~T7ll : 

[H(x),[H(x),A 1 ]] + 3[H(x),A 2 ]=0, 

[H(x),[H(x),A 2 ]] = 0. (3.38) 

Proof. Despite it might be tempting to attack the problem similarly as in the case of 
Lemma 1, a direct calculation reveals that one cannot avoid binomial expressions with 
double summation involving linear combinations of quadratic terms (products of two 
binomial coefficients), which are extraordinary difficult to deal with. Fortunately, as 
we demonstrate below, there exists an elegant algebraic recursive procedure originating 
from an extra symmetry of the generator H(x). Since we are dealing with quadratic 
expressions in H(a;), whose blocks Yi^ a \x) are singular with one- dimensional null-space, 
additional information about null-vectors of (H* )) 2 will be required as well. 

Here we shall label quantum space components with the Pauli basis. According to 
the structure (13.9113. 1Q|I3. lip the Eqs. ( 13. 38ft are equivalent to five identities which can 
be cast in terms of operators over "H a ® "H a , {D^, D^ }, s G {0, +, — , z} (temporarily 
dropping spectral parameter dependence for the rest of this proof), 

D? :=[H,[H,A?]] = 0, 
Bl :=[H,[H,A z 1 ]] = 0, 
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D+:=[H,[H,A+]]+3[H,A+] = 0, 

D7:=[H,[H,Ar]] = 0, 

D+:=[H,[H,A+]]=0. (3.39) 

Consistently with our previous notation, we will place additional subscript index a, e.g. 
Ag , when referring to a single a-subspace. 

The key ingredient here is notification of a 'conserved charge' AT/, 

[H,A-] = 0, (3.40) 

connecting two adjacent a-blocks, 

H (a+1) A W- =A W- H (a)_ ( 341 ) 

Because the above identity is to hold regardless of the value of x, we essentially have to 
prove [X( Q ) p , A" ] = for all residue matrices (1A.5I) for p — 0, 1, ... a, demanding to 
verify the identity 

(I + 1)X&? - kxi% - (a - J + 1)X^ 1)P + (a-k + l)x£f = 0, (3.42) 

for every k = 0,1, . . . ,a + l, I = 0,1, . . .a. In fact, it is sufficient to consider the identity 
expressed in terms of tensors Y^ p by virtue of parity symmetry (1A.8I) . 

(I + l)<m )P - kY«% -(a-l + l)r fc ( J +1)p + (a - k + 1)Y^ P = 0, (3.43) 

which reduces to trivially verifiable combinatorial identities upon substitution (1 A. 71) . 
Next we state algebraic relations among A-[ a 2 , which are straightforwardly verified 



using explicit representations ( 13.30fl3~3"4"l) . namely 

A (a+D0 A (a)- = A (a)- A W0 < 

A (a + l) ZA (a)- = A («)- A (a)« _ 2A (a)- ; 

A (a) + A (a)- = A (o-l)- A (o-l)+ + A («)z 5 

A ( a) + A ( Q) _ = A (-l)- A (-D+, (344) 

The idea is then to derive recursive relation in a for the operators {D^ , D!> }, 
and use induction in a, along with the trivial initial conditions D^ a 2 = 0, for a = 0,1, 
which are easy to check (e.g. by direct evaluation), to prove the identities 1 13. 391) . Since 
all obtained recursions are treated in analogous way, we choose to work out explicitly 
the one with D^\ After expanding double commutator, 

D (o.)+ = ( H («)) 2 A^ )+ - 2H^A^ )+ H( a+1 ) + A^ )+ (H( a+1 >) 2 , (3.45) 

multiplying by A^ from the right, using (i) H^ a+1 ^A^ = A[* H^ and (ii) 
Ag A-i — ^i -^-2 > an d commuting A{* to the left, we obtain 

D («)+ A («)- = a$*- 1) ~T>£- 1)+ . (3.46) 



The relation we have just established is however not enough to conclude on vanishing of 
Df> , provided Dg = 0. The reason lies in the non-invertibility of the rectangular 

matrix A[* . We can nonetheless cure this weakness if we show that there exists an 
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additional (a + 2)-dimensional vector u( a+1 ), linearly independent of the column space 
of A® , which is in a null-space of D;, , 

D («)+ u («+i) = o. (3.47) 

Notice that the remaining a + 1 columns of A{* are indeed linearly independent 
(another simple calculation). 

A crucial observation is, that for every a-sector, there exist a unique pair of (a + 1)- 
dimensional null- vectors u^.v^, 



v< a > 



(i,-i,i,-i,...) T = E(- 1 )*i*>» ( 3 - 48 ) 

k=0 

a 

u<«> = (0, -1, 2, -3, ... , (-lTaf = ]T(-l) fc A; \k) , (3.49) 



fc=0 



such that 



H (a) v (a) = o, (HW) 2 u( Q ' = 0, H( Q >u( a » = -vW (3.50) 



a 
x 

In order to prove our case fl3.45|3.47l) . we require another set of identities, expressing the 
action of a squared generator a-blocks (H^-*) 2 on the vector vS a \ and transformation 
of both u^ a ) and v^ a ' under the action of A^ 2 . For the sake of brevity, we entirely omit 



their justification here (it can be found in Appendix C): 
A («)o v («) = A (")+ v («+D = o, 

A< o) V a > = -2ocv&\ 

A W+ v («+i) = av to (3.51) 

and 

A W+ u (a+l) = 2v (a)_ (3 _ 52) 

Accounting for expansion of D 2 (13. 45 p . and using auxiliary identities (I3.5ip .f l3~52|) . we 
show D 2 u < ' Q+1 ) = 0. Linear independence of u(° +1 ) from the column space of A" 
therefore allows for its extension to an invertible matrix A t a by adding u^ a+1 ^ as the 
(a + 2)-th column vector, yielding 

D ( Q)+ = Aj^-D^^A^)- 1 , (3.53) 

from where it immediately follows D 2 = if D 2 = 0- 

An entirely analogous reasoning applies to the remaining four cases from ( I3.39p . 
Using null-vector of (H^) 2 we derive the action of the operators {A* a ] for s = {0, +, z}, 

A («)o u (a) = _ 2av («) j (3.54) 

A («)* u (a) = _ 2av H - 2(a - 2)u( Q \ (3.55) 

which justifies adding u^- 1 (or u^ a+1 ^ in the case of D^ ) to the columns of Aj* 
when operating by the corresponding D^ from the left. Essentially it sufficient to 
demonstrate that D^ preserve the null-space ker(H 2 ). 
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In the case of diagonal blocks D^ and D^ , after multiplying them by the 
conserved charge A" from the right, we use identities (13. 44ft and 

H («) A (-i)- = A («-i)- H («-i) ( 3 _ 56 ) 

to bring A" to the left, finishing with 

Bl a \[ a) ~ = A^D^, a e {0, z}. (3.57) 

An extra linear term in D^ is of no importance, as it cancels out regardless of its 
prefactor. Finally we take care of D^ , arriving at the following coupled operator 
recursion 

D («)+ A H- = aS q - 1) -dS q)+ + Di a)z . (3.58) 

We have nevertheless already proven that Df* = for every a, hence the recurrence 
becomes of the same type as the ones above. 

4. Properities of the exterior integrability structures 

4-1. Properties of the R-matrix 

i?-matrices are required to obey additional compatibility-type condition (the celebrated 
Yang-Baxter equation, or in our notation, the braid group relation) imposed on a triple- 
product of auxiliary spaces *Hf 3 , 

(l®R(A,//))(R(A,?7)®l)(l®RGu,77)) = (R{n,r])®t){t®R{\,r])){R{\, (j,)®t), (4.1) 

which automatically ensures associativity of intertwining property over multiple spaces 
% a . Despite in this paper we only strictly prove the RLL relation ( I2.13P we discuss in 
Section [5] some other related results [8j [H] from which (14. ip should also follow. 

Unlike in most often encountered cases of integrable models (e.g. in fundamental 
models), the .R-matrix here is not of difference type, i.e. its elements do not depend on 
the difference of the involved spectral parameters only. Yet, the difference of spectral 
parameters, curiously enough, enters in a way ( 13. ip which is reminiscent of a Lie group 
structure. 

Additionally, one observes the following interesting properties of the R operator (all 
following directly from explicit representation R(A, fi) = exp ((A — /x)H ((A + ju)/2)) and 
the properties of the generator H(x)). 

(i) Regularity: 

R(A,A) = 1. (4.2) 

(ii) P -symmetry [See definition (I3.19p ]: 

7r a (R(A,/i)) = R(/i,A). (4.3) 

(iii) Orthogonality: 

R(A,aOR(>,A) = 1. (4.4) 
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(iv) From (ii) and (iii) we immediately derive another nice property, namely that the 
eigenvalues of the matrix PR(A,/i) can only be ±1, namely 

(PR(A,^)) 2 = 1 (4.5) 

(v) All eigenvalues of R(A, fx) are equal to 1. The R operator has therefore a non-trivial 
Jordan canonical form, where each (a + 1) x (a + 1) matrix R( a )(A, /i) is similar to 
a single irreducible Jordan block [following from (1A.12J) . 

(vi) The matrix R(x + \y,x — |) is holomorphic in both x,y G C except at x G |Z + 
where it has simple poles. This follows from the fact that the generator H(x) has 
simple poles at x G |Z + as well, property (v) which terminates the exponential 
series after a + 1 terms in subspace nC\ and a curious nilpotent algebra among 
its residua X^ = Res 2 , = pH(x) (1A.7[) . namely 

X (p) X (m) =0 if p > m. (4.6) 

The property (14.61) can be studied in each space %& separately, where it can be 
proven by application of inductive arguments on elementary binomial identities. 

(vii) At the poles, actually, where A + /i6 Z + , the intertwining of the product of Lax 
operators can be implemented by taking a residuum of the RLL relation (12.131) . 

(viii) Comparing the representation (I2.6P of sl(2) with the transposed onqjj], which should 
be equivalent 

(-l) s Aj(A)=U(A)A_ s (A)U- 1 (A), sG {+,0,-}, (4.7) 

where U(A) is a diagonal operator from End("H a ) (invertible for A $■ \l> + ), 

U?(\) = 5 k ^y (4.8) 

we obtain the corresponding transposal symmetry for the i?-matrix 

(U(A) ® U(aO)R(A, /iXU-V) ® U^ X (A)) = R T (/i, A). (4.9) 

Note that this is a kind of Liouvillian PT symmetry of the type proposed in Ref. |25j . 
The sign factor (— l) s in ( 14. 7p is a consequence of non- canonical (real) representation 
of flMD- 

4-2. Properties of the monodromy matrix 

Rich structure and properties of the R operator discussed above are also inherited by 
the corresponding exterior monodromy operator T(A) or its matrix elements (j2.22p . We 
list some of the most remarkable properties that we have observed here, with the hope 
that they will find useful future applications (e.g. those discussed in Section [5]). 

Firstly, for a given system size n, selection rule (I2.24p implies that the monodromy 
matrix is banded, i.e. 

T l k (X) = 0, if \k-l\>n. (4.10) 

| Transposition is defined, as usual, (|fc) (l\) T = \l) (k\, without complex conjugation. 
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Furthermore, we claim that only the elements T z fc (A) from a (n + 1) x (n + 1) square, 
namely for k < n, I < n are linearly dependent physical operators. Therefore, for a 
fixed distance from the diagonal q = \k — l\, only n — q matrix elements are linearly 
independent, while all others can be expressed in terms of those 

n— q n—q 

fc=0 fc=o 

where c^ qlk {\) are some rational functions of A with integer coefficients. 

Secondly, we were looking for linear combinations of magnetization (particle- 
number) preserving diagonal matrix elements T"*(A) that would form a commuting 
family. Up to linear dependences we conjecture (based on empirical evidence) that 
there exists a single commuting linear combination besides 5(A) = T® (A), namely 

*w = B-ir + '(:)^^(A), (4.12, 

so that 

[5(A), 50)] = 0, [5(A), 5^)] = 0. (4.13) 

Thirdly, as in our problem we are dealing with non-normal operators one may want 
to understand the connection between the transposed ([Tl(X)} T ) and original (T/^A)) 
monodromy elements. Writing the reflection parity operator in quantum spin space 
Q = Q^ 1 E End("Hf n ), Q \v 1: v 2 , ■ ■ ■ , u n ) = \u n , v n -i, ■ ■ ■ , ^i), we find immediately 
[applying Eqs. A223KT} ' 

K(A)) T = (-1)^ ^ ^QTfWQ, (4.14) 

or more compactl}l5j, writing a partial transpose with respect to %f n by superscript T s , 

T Ts (A) = U(A)gT(A)gU- 1 (A), (4.15) 

were U(A) = diag(l, — 1, 1, — 1 . . .)U(A). Furthermore, the reflected monodromy 
elements QT l k (X)Q can be in turn expressed in terms of linear combination of T£_ k , 7 -(A). 
For example, we state the connection explicitly for 00 matrix element 

\ / k= Q \ / 

Finally, let us consider the action of monodromy elements on spin states with a 
small number of quasi-particle excitations. Let |O m ) G 7if n denote an arbitrary state 
from ( n J dimensional subspace with exactly m spins-up (and all other n — m spins 
down), i.e. M \fi m ) = (2m — n) \Q m ), and let \Cl m ) denote a state with m down-spins, 

% Yet, another form of PT-like symmetry [25] , 
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\fi m ) = \Vt n _ m ). Direct inspection using explicit representation (I2.8P reveals the action 
on the vacuum state in terms of a shift of spectral parameter 

T/ +g (A) |n ) = ( 2A g " Z ) ( 2A ~ 2/ ) V°(A - k) |fi > , (4.17) 

T/ +9 (A) |Oo) = (' + 9 )t*(A - fc) |O ) . (4.18) 

whereas, one can write similar but more general expressions for the m-particle sectors 

T l l+q {x) 1^) = ^Tr«f (A)T fc fc +g (A - (/ - k))) |n m ) , (4.19) 



fc=0 
m 



r/+V) |n m ) = £ s^{\)T k k+q (a - (/ - fc)) |n m ) , (4.20) 

fc=0 

where q E Z + can be interpreted as the number of quasi-particles created and I > m 
for the relations to be non-trivial. Remarkably, the rational functions rf™{X), sf™(\), 
again having integer coefficients, and only simple poles at |Z + , do not depend on system 
size n. For the purpose of treating the NESS density operator, say for developing an 
ABA procedure for diagonalizing it, it should be handy to control transposed matrix 
elements at negative spectral parameter —A, corresponding to Hermitian conjugation at 
real value of dissipation s = 2i/A, for which A = —A and Sft(\) = S T (—X). Let us write 

f l k (x) := (-ir[T*(-z)f, (4.21) 

where a sign factor (— l) n is put for convenience. Then, straightforward inspection again 
suggests remarkable connections: 

in 
T l l+q {\) |n m ) = £ f?;rWT k k +q (A - (S + I + k)) |n m ) , (4.22) 

fc=0 

III 

f/ +9 (A) |n m ) = Y,9f;rWT k k +q (a + (q + 1 + k)) |n m ) , (4.23) 

k=0 

III 

T^W |n m ) = $> 9 f (-A)f fc fe +(? (A - (g + I + k)) |n m ) , (4.24) 

fc=0 

III 

f l l+q {\) |n m ) = £ /£ m (-A)T fc fe+9 (A + (q + I + k)) |n m ) , (4.25) 

fc=0 

where relations are already non-trivial for any /, and integer coefficient rational functions 
// 9 fe m (A) 5 sf)T(A) again do not depend on size n. 

5. Discussion 

After this work has been completed, we have learned about Refs. [U [HI [7] where 
related infinitely-dimensional i?-matrices have been constructed using manifestly sl(2)- 
symmetric form of Lax and i?-matrices. It seems that such a universal sl(2) i?-matrix 
might be useful in the context of QCD and high-energy physics whereas in condensed 
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matter physics the non-unitarity of the general infinitely-dimensional representation 
seems to be only compatible with phenomena far from equlibrium which we discuss 
here. 

In fact, our Lax matrix f 1 2 . 1 j) becomes s[(2)-symmetric after multiplying by a z , 
L(A) = L(A)cx z = B+(A) (8) a' + B_(A) <g> a + + B Z (A) <g> cx z = B(A) ■ a, where 
B + = A_,B_ = — A + ,B Z = A are canonical generators of infinitely dimensional 
representation of sl(2) with representation parameter A Q The i?-matrices resulting 
from solving RLL relations for the two forms of L-matrices, L(A) and L(A), are different 
but closely related. Nevertheless, the results presented in this paper are more explicit 
and detailed in connection to a different form of a transfer matrix as they are taylored 
for non-equilibrium condensed matter applications, and hence they are essentially non- 
overlapping with those of Ref. [Hj. Although the sl(2)-symmetric L-matrix generates 
a related transfer matrix, namely S(\)(a z )® n , and yields an identical NESS density 
operator S(X)S^(X), we have a good reason to use also a symmetry-broken representation 
of the Lax matrix. Namely, only in our representation the MPA for S(X) generates a 
convergent sum of local operators [13] in the g-deformed case of anisotropic XX Z model 
(see discussion below, in subsect. 15. 2p . 

We foresee two immediate interesting applications of the exterior (non-equilibrium) 
integrability formulated here. 

5.1. Algebraic Bethe Ansatz and spectrum of the density operator 

A tempting proposal following from our construction is the construction of ABA 
procedure for diagonalizing NESS density operator. This could be particularly 
interesting in the light of recent suggestions [201 I2S] that the spectral properties 
of equilibrium and non-equilibrium density operators can be used as indicators of 
integrability (or exact solvability) similar as in the idea of quantum chaos. 

Algebraic form of Bethe ansatz allows for construction of an eigensystem for a 
family of mutually commuting transfer operators. The procedure is based on the quasi- 
particle modes created under the action of (off-diagonal) elements of the monodromy 
matrix T(A) = L(A)® sn . Many-particle excitations arise as a string of monodromy 
elements, operating on a specially chosen reference state. The role of i?-matrix is to 
prescribe quadratic algebraic relations among elements with different value of spectral 
parameter (which are interpreted as quasiparticle momenta), enabling for construction 
of eigenstates of the quantum transfer operator. A set of n spectral parameters {A^} for 
n-particle excitations has to be chosen accordingly in order to eliminate unwanted terms 
(those that are not the eigenvectors) which unavoidably emerge during commutation of 
the elements of T. The latter condition gives rise to famous Bethe ansatz equations [I]. 

As construction of ABA in this case, due to Cholesky structure of the diagonalizing 
operator H2.4J) . does not seem to be straightforward, we outline here only the first 
step. Namely on how to obtain single quasi-particle excitations, i.e. eigenvalues and 

+ In the usual complex representation we in addition have to re-define the generators Bj- — > — iB±. 
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eigenvectors of p^X) = S(\)S T (-\) = (-l) n T °(A)f °(A), where A = 2i/e e iM of the 
type |fii). Applying the connections fl4.2204.23p and the RTT relation f l2TT5|) for Hi a=1) 
sector only, therefore using only the 2x2 block R^, we arrive at the useful identity 

(-l)"j5oo(A):Z?(ji) |flo> = [t(A)] 2 A(A,ji)2?(/*) |fi > (5.1) 

//((A + /i - l)t(A)t(/x) - 2A(p - A)t(A + !)*(// - 1)) 

+ (,-x)(x-, + i) Tl (A) |fio) 

2 / iA(A + l/2)t(A)t( / i-l) 
+ (A + lKA-,, + 1) T ^ A + 1 )I^' 
where t(A) := A n and A(/i, A) is a quasi-particle dispersion relation: 

(A + ^)(A + /i-l) 

A/^,A :=— — — — . (5.2) 

{\-fi){X-fi + l) 

There are two single quasi-particle states with the same eigenvalue A(//i, A) = A(/i 2 , A), 
which can be parametrized in terms of a single complex rapidity parameter £, as 
//i = o(l + (A + 1)0; ^2 = |(1 + (A — l)/£). Hence the single-particle ABA is already 
a nontrivial combination of two terms 

I*) = (CilW) + C 2 7?Gu 2 )) |O ) (5.3) 

As the three operators on the RHS of (15. ip are linearly independent, the requirement 
that the two unwanted terms, proportional to vectors T°(A) \flo) and T°(A + 1) |f2o), 
cancel, i.e. to have Poo(A) |\l/) = A |\1/), results in requiring that a 2 x 2 system of 
equations for C\,Ci has a nontrivial solution, i.e. 

!-(A + l)e W e + A-l V / l-A (A + l)g + A-l 

i + (a + i)J u-a+v vi+d (A+i)e-A+r l - J 

This can be understood as a Bethe equation for single particle eigenvectors of NESS, 
with eigenvalue A(~(l + (A + 1)£), A). 

However, generalizing this procedure to multiple excitations seems far from trivial 
and should be a challenge for future work. 

5.2. The anisotropic XX Z model and a new family of quasi-local conservation laws 

As has been pointed out in Ref. [H], the infinitely dimensional i?-matrix also exists for 
continuous representations of the quantum group U q (si(2)) hence all our constructions 
of exterior integrability should be q— deformable and should translate to the boundary- 
driven anisotropic XX Z spin chain [221 E31 [15] where the Hamiltonian density h in 
(12. 2 p should be replaced by h = 2a + ® a~ + 2cr - <g> a + + Acr z <g> a z with the anisotropy 
parameter A. 

Most interesting there is the question, whether the recently discovered quasi- 
local conservation law [22] can be generalized and extended to a whole family. 
In integrable theories, local conserved quantities are usually obtained in terms of 
logarithmic derivatives of transfer matrices around some trivial values of the spectral 
parameter. Here, the spectral parameter is non-standard and is related to coupling 
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to the environment, hence the derived conserved quantities can have different spin-flip 
symmetry K = (<j x )® n as in the standard case [9J [16] where due to the equivalence 
of quantum spin and auxilliary spaces the K symmetry is imposed to the L- and R- 
matrices as well and henceforth to all so-derived families of conservation laws. In the 
exterior integrability problem, however, the A-symmetry is explicitly broken, resulting 
in (potentially quasi-local) conserved quantities which may yield non-trivial Drude- 
weight bounds [13] even in the absence of external magnetic field. 

For example, writing for the moment the commuting transfer matrix 1 12.llj2.12p as 
a function of dissipation e (in notation of Ref. [23J), which is polynomial for finite n, the 
conservation law Z which has been proposed and implemented in Ref. [22] and which is 
quasi-local for |A| < 1 is simply Z = ^S(e)\ E= o. We conjecture that a further tower of 
(quasi-local, in case |A| < 1) conservation laws which break the A-symmetry is given 
by higher logarithmic derivatives 

Z * = ^=i lo gS(£)| £ =o, k = 1,2,... (5.5) 

The even order logarithmic derivatives vanish as a consequence of an interesting identity 

S-\e) = S(-e), (5.6) 

which can be easily proven. Details on these constructions shall be presented elsewhere. 

5. 3. Conclusion 

We have provided a new link between the matrix product ansatz and Yang-Baxter 
integrability in the context of non-equilibrium quantum physics, which is fundamentally 
different than the one which exists on the level of closed quantum systems pp. The 
first fundamental difference is in the role of spectral parameter of the integrable 
theory which is now taken by a continuous representation parameter of infinite- 
dimensional representation of the underlying quantum symmetry of the model. The 
second fundamental difference is the formulation of the transfer matrix, which is here, 
due to infinite-dimensionality of the representation space, taken by the ground-state 
expectation instead of a trace. Generalizations to other quantum integrable models seem 
straightforward, the most obvious one being perhaps the multi-component quantum 
hopping model [3~T] . 
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Appendix A. Explicit expression of the generator H(x) 

Here we provide three additional, trivially equivalent, but useful forms for the a-block 
of the generator H^ a \x) of the exterior i?-matrix. 

Compact form. 

and (13 .3p elsewhere, where use the C-number extension of the binomial symbol 
'x\ x(x — l)(x — 2) • • • {x — k + 1) 



,*;■ fc(ife-i)(jfe-2)...i ' xeCm (A,3) 

Showing that (1A.2J) is identical to the expressions (13. 3p amounts to checking that they 
have identical residua. This in turn leads us to the third, 

Residue form. For most purposes we find useful the following matrix-valued residue 
decomposition of H^, with poles positioned at x = m/2, m = 0, 1, 2, . . .: 

a 

H^(x) = J2x (a)m fm(x), (A.4) 

m=0 

X (a)m :=Res x=f H^(x), (A.5) 

and matrix coefficients X^ m are given in P-symmetric form explicitly as [see (13. 3p ] 

Y (a)m _ f ( v (a)m v {a)m \ /a m 

r(a)m __ ( 1 Nfc-m-l/ ™ \ / A ' l x 



where 9 X := 1 if x > and 6^ := if x < 0. Note that in terms of parity operator H 3 . 2 1 f) : 
X = -(Y-PYP). (A.8) 

Jordan form. Here is the form which is in fact equivalent to a Jordan decomposition 
of H^. Let W a (x) be an upper triangular (a + 1) x (a + 1) matrix with entries: 

which vanish if fe > Z, and A a a strictly lower triangular matrix with constant entries: 

A h = j—T if fc>/ ' (A ' 10) 

A ti = ° if fc < L ( A - n ) 

Then, we have the following decomposition: 

H a (x) = W Q (x)A a (W Q (x)) _1 , (A.12) 

where, noting, (W a (x))~ is again upper triangular. 
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Appendix B. Verification of the HLL relation. 

By utilizing the residue decomposition (I A. 5 1) of a-block H^ a \x) in terms of X( a ', we 
calculate component-wise expansions of 

[HW(i),Ao(i)]=Ai (B.l) 

using definition ( 13. lip with explicit form of MPA amplitudes (12. 8p . Let us initially 
consider the diagonal physical components E 00 (recall that E u is automatically obeyed 
by virtue of symmetry (13. 28ft ). where block operators Aq preserve hC . It is also 
noteworthy that the corresponding A" on the right side is bidiagonal (constant) 
matrix, therefore a direct calculation leads to, after isolating matrix coefficients in front 
of every simple pole f p (x), p = 0, 1, . . . , a and for every matrix element k, I = 0, 1, . . . , a 

(k-l)(k + l- a)X$ p + ( a -l + l)(l-l - P )xi% 

+ (k-a)(k-p)xi% = 0. (B.2) 

When off-diagonal physical components are considered, two adjacent sectors a and 
a + 1 will get coupled, i.e. for arbitrary a we evaluate 

[H< a )(x), A { a)0 \x)} = A[ a)0 \ [B»(x), A { a)w ] = A[ a)w , 

or in component notation, 

(I - p - 1)(2I +p - 2a - 2)xi% + (2l-p)(l+p- a)X^ 

_(2fc-p)(fc + p -«)xg +1) ' p -(fc-p)(2A; + p-2a)A£t! ) ' P = 0, (B.3) 

k(2k + p-2a- 2)XI% + (2k-p)(k-a- l)xg )p 

-(21 -p)(l-a- 1)X£ +1) ' P -(l + l)(2l + p- 2a)X^ = 0, (B.4) 

for physical components E 01 and E 10 , in respective order. All the relations above are 
of homogeneous kind because the singularities {f p (x)} are absent in the right-hand side 
of the RLL relation (involving Ax). It is worth noticing that parameter p enters in the 
amplitudes from 'fusion' with poles f p (x), by virtue of partial fraction expansion, 

w(a,k,l) -2x _ w(a,k,l) -p _ 
x — p/2 x — p/2 

For instance, working out ( 1B.2I) explicitly at fixed a and for the pole f p (x) (for sake of 
clarity we omit x-dependence from the amplitudes), we have 

{[Xl°*f p , K\x)}) kl = f p \(ay Q ^ - aU4)xif + aUa^X^ - a^_ l0 +X^ 



k+l,l 



= f p (k -l)(k + l- a)X { k f + f p (l-l-p)(a-l + l)xi% 

+ f p (k - a)(k - P )X { k % 

-2[(a-l-l)xi% + (k-a)xi%], 

Beside matrix-residue part, non-singular terms are produced as well. By collecting 
together non-singular contributions from all {f p {x)} and matching them to non- 



Exterior integrability 23 

vanishing (linear in x) terms on the right, we get an additional set of conditions which 
are to be satisfied: 

a— 1 1 

p=0 

at the E 00 component, and analogously 

E {4% + 4? - Xg+ 1)p - X<S$>} = 0, (B.7) 

p=0 



E { ( 2a " l ~ P + ^ X i% + V-P + *)4T +(p~k- a)X^ 1)p 

p=0 

+ (k + p- 2a)X^+f } = (3k - a)5 Kl + (3A; - 2a)4+i,*, (B.8) 



E { kX i% "(*-«- m { k f + (l~a- 1)X£+ 1)P - (I + l)X^ p } 

p=0 

= {a-l + l)6 k>l -(l + l)6 kil+1 , (B.9) 

at E 01 (equations f ]B.7j) .( )R8|) ) and E 10 (equation (1B.9J) ). However, as the latter set 
of expressions is rather tedious for further analytical manipulations, we decide at this 
point to take a different (however equivalent) strategy and rather employ the first form 
of the generator (I3.3P . 

For the sake of compactness, we shall only provide explicit calculation to justify 
validity for the set of equations pertaining to non-singular part for physical components 
E 00 and E 01 , whereas an entirely equivalent procedure applies to show the identity 
associated with E 10 component. We start with the diagonal element, where from fIB.lj) 
we obtain 

( \hM(x), A[ q)00 1) = (k - l)(k + I - a)H^(x) 

V L j / k,l 

+ {a-l + l)(l-l-2x)H i k "l 1 (B.10) 

+ (k-a)(k-2x)Hl% 

= {a- 2l)S ktl + 2(a-l+ l)5 k+lil = (a^ 00 ) . 

\ / k,l 

We introduce diagonal index 5 := k — I and focus initially on situation 8 > 1, where 
equations become homogeneous. Projecting out components coupled to any f p (x) (as 
they are irrelevant for this part) we find the requirement 

m=0 v / \ / 

'l + 5+l\ f 8 



n + S-a) =0, <5>1, (B.ll) 

/ / \m. 
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which is obviously true for all/ = 0, 1, . . . , a, based on a well-known binomial identity 



.m, 

m=0 N ' 

The diagonal cases follow after plugging 5 = (beware of all the corresponding 
prefactors), where only m = contributes, yielding 

l [ \a - I + 1) - ( l + X Vz - a) = a - 21, (B.13) 

which correctly reproduces diagonal elements of the right-hand side of (IB.lOp . The same 
argument of course applies when S := I — k > 2, the only difference being the indices in 
the poles get reversed, i.e. the amplitudes fuse with the element H^_ k a _ l . This however 
leads to the same argument based on the identity ( IB. 121) as long as the non-singular 
part of the expression is considered only. Thus it remains to be checked in the case 
when k = I — 1, where 



r(«) 



(2l- a -l)H^> a ^(x) + (a-l+m-l-2x) H£>(x)-H} 



r( a ) („\ Tj( a ), 



X 



= l(2l-a-l)(a-l + l)f a -i - i( a - / + 1)(/ - 1 - 2x)(/,_ 1 + / a _,) 

= 2(a-/ + l), (B.14) 

where again the correct result 2(a — I — 1)^+1,/ is reproduced. 

The same procedure applies for the off-diagonal physical component, where element 
from two neighboring a-subspaces are involved - e.g. for the E 01 we have to show that 

(Z - 1 - a + x)(l - 1 - 2x)H { k a l 1 (x) + (a - I - 2x)(-l + x)H${x) 

+ {a-k-2x){k- x)H { k a l +1 \x) + (a - k - x){k - 2x)H^f 

= (3k - a)8 kjl + (3k - 2a)8 k+1>l . (B.15) 

By focusing once more on a non-singular part after resolving expansion in terms of 
{fm( x )} an d terms which are now linear functions in x, we find the vanishing of the 
latter is implied, for S = k — I > 1 and for I — k > 2 based according to (IB. 12|) in 
conjunction with another identity 

£ ( - irm Cl) =o ' (Ri6) 

m=0 ^ ' 

regardless of the form of corresponding prefactors (which are functions of parameters 
a, k, I). It is left to check for special cases now - at k = I we calculate 

(I - a - 1 - x)(l - 1 - 2x)H<if_ l (x) + (a - I - x)(l - 2x)H^f(x) 
+ (a-l-2x)(l- x) (< +1) (x) - Hif(x)) 

= - 1(1 -l-a + x)(l-l- 2x)fi-i + (l + l)(a-l-x)(l- 2x)fi 
-(a- I- 2x)(l - x)f a -i] =31 -a, (B.17) 
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and for k = I — 1, 

(a-l + l-x)(l-l-2x) fer +1) (a:) - 

+ (a-l- 2s)(-J + x)H^\ }l (x) 

+ (a-l + l-2x)(l-l- x)H^u\x) 
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#. 



(a)(x) ' 



(«-/ + l-x)(/-l -2x)/j_i 
- (a - / + l)(a - / - 2x)(-/ + x)f a -i 
-(a-l + 2)(a - / + 1 - 2x)(Z - 1 - x)/ a _ z+ i 
= 3(7-1) -2a. 



(B.1S 



One applies the same arguments to show the remaining case of the E component. 

Appendix C. Nullspace vectors of H(x) and H(x) 2 

1 . Vector v( a ) is in the kernel ofH^ a \x). It is sufficient to prove that 

X (a)p v (a) =0? p = 0,l,...a. (C.l) 

This requirement is in fact implied by two separate (stronger, i.e. sufficient) conditions 



Y(o)p v («) 



» 



p v («) 



a v ( a ) 



•1 Q v 



(C.2) 



The second being obviously satisfied, we focus on the first one and employ the component 



notation v 



(a) 

i 



E 

1=0 



■l) 1 . We have to show that 

'k\ fk-l-l 

l)\ p-l 



\k— p— 1 



vfc-1 



which can be in turn recast into 



Em: 



P+i 



z=o 



k\ fk-l-l 



p — I 



(C.3) 



(C.4) 



One can quickly check that for p = the contribution comes only from / = and the 
equation trivially holds, which serves as our basis of induction. Next, by induction step 
we move to p — ¥ p + 1, which after application of the Pascal's rule yields 

p+i 



B-i: 



p+z 



£(-W 



v 



k\ fk-l-l 

p — I 



\p+i 



k-l 
p-l + 1 



(C.5) 



'jfc - / - 1\ / fc - / 

P-l )~ \p - I + 1 , 

p+i 

On the right-hand side we retrieved an expression from the previous step plus an extra 
sum. Introducing a summand function V(7, 1), it is thus necessary to show that 

/jb\ a - r 

,7-1 



Enr.o ==»-!)' (J) (! 

Z=0 2=0 \ / \ 



0, 7>1, 



(C.6) 
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We rely on the observation that the sum at hand is Gosper- summable [2TJ, i.e. because 
the summand obeys the following recursive formula 



V( 1 ,l) = A l 



— )n-./) 

7 



7^0 



(C.7) 



using the definition of the forward difference operator A m A(m) := A(m + 1) — A(m), 
the resulting telescoping series with finite support vanishes. 

2. Vector u^ = (0, -1, 2, -3, ... , (-l) a a) T is (i) an eigenvector ofY^ p and tt^Y^p) 
with eigenvalue —1, provided p > 1. Additionally, (ii) for initial value p = we have 
Y(«)° u ( a ) = and 7r a (Y( Q )°)u( a ) = -av^. Therefore, (i) and (ii), together with 1. 
imply that (H( a )) 2 u( Q ) = 0. 

Initially for p > 1, rewriting the action of Y^ p in components, and accounting for 



u| = (—1)'/, we obtain 



J2y^ a) = J2(-i) k - p - l ( 

1=0 1=0 ^ 



xfk\ (k-l-l 
l){ p-l 



■i)H 



■i) k k. 



(C. 



Since for k = it is evidently valid, we divide by k and by means of (, ,) = !(,) and 
reformulate it as 

z=o i=o v / \ P / 

Beginning with the basis of induction at p = 1 , we first show 

i 

5>F(p,Z) = -F(l,l) = l. 
z=o 

Proceeding with the inductive step p — > p + 1 we find 

p+i p+i 



(CIO) 



k-l-l 
p-l + l 



£F (P+ u)=£(-i)<-*+'r: )( 

p+1 

= £(-!)'- 

P /'"i-i- /r i 

=E^.')-E(- i ) , -' , n_i 

«=o z=o ^ 



k-l 

/-l 

p+1 



k-l-l 

p — I 



k-l 
p-l + l 



k-l 
p-l + l 



It is necessary to show, that the second sum always vanishes for 7 > 2, 

Z=0 i=0 v / \ / / 

which is again summed up by help of the recursive formula for the summand, 



Vfr,l) = A l 



LI-7 



(C.ll) 



(C.12) 



(C.13) 
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Finally, we prove exceptional cases at p = 0, where only I = contributes, and 
consequently 

(Y^°n^) k = Y$°u^ = 0, 

(*.(¥(«>><«>), = y£i>W = (-l) fc (-«) = -c4 Q) . (C.14) 

Hence, 7r a (Y < - Q ^ )0 )u^ = -av'"' and consequently also 

X H0 U («) = « (a)_ (Q 15 ) 

By combining the above results with the residue form of H^ a '(x) (I A. 71) we conclude 

H< a '(x)u (a) = -v (a) . (C.16) 

x 
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